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Abstract. We study the statistics area, bounce and dinv on the set of parallelogram polyomi- 
noes having a rectangular m times n bounding box. We show that the bi-statistics (area, bounce) 
and (area, dinv) give rise to the same q, i-analogue of Narayana numbers which was introduced 
by two of the authors in 1 3j. We prove the main conjectures of that paper: the q, t-Narayana 
polynomials are symmetric in both q and t, and m and n. This is accomplished by providing a 
symmetric functions interpretation of the q, t-Narayana polynomials which relates them to the 
fT") famous diagonal harmonics. 

o 

1. Introduction 

03 

A parallelogram polyomino having an m x n bounding box is a polyomino in a rectangle 
consisting of m x n cells that is formed by cutting out two (possibly empty) non-touching 
standard Young tableaux which have corners at (0, n) and (to, 0). An example of a parallelogram 
polyomino having a 12 x 7 bounding box is illustrated in Figure [Tj Let Polyo mn be the set of 
\mJ all parallelogram polyominoes having a rectangular to x n bounding box. The cardinality of 

u Polyo m n is known to be N(m + n — 1, m) where for positive integers a and b, 

i— i are the famous Narayana numbers. Two authors of this work introduced [3] two statistics on 

these combinatorial objects, area and bounce, which led to a q, i-analogue of the Narayana 
j>. numbers N(m + n — 1, to), namely 

g Nara m , n (g,t):= £ g area(P) f bounce(F) ) 

OO PePolyo min 

that they called the q, i-Narayana polynomial. In that same work it was conjectured that the 
q, t-Narayana polynomials are symmetric in q and t, and as expressions were also symmetric in 
to and n. We introduce a new statistic dinv which gives a new q, t-analogue of the same numbers 



(N 



X 



Nara ro ,„( 9 ,t) := ^^inv(P )f area(P) 

p 



The following theorem establishes a relation between these two polynomials. 



Theorem 4.1 , For all m > 1 and n > 1, we have 

Nara TO) „(g,f) = Nara„ im (g, t). 

We give two proofs of this result, one using an explicit bijection, and another one using a 
recursion. The main result of this paper is the proof of the symmetries conjectured in [3]. 



Theorem 6.2, For all m > 1 and n > 1, we have 

Nara mi „(g,t) = Nara mjn (t, q) 



Key words and phrases, q, f-Narayana, parallelogram polyominoes, parking functions. 
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f Supported by NSF grant DMS=0800273. 
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and 

In particular 



Nara mj „(g,t) = Nara n , m (g,t) 



Nara m ,„(g,f) = Nara mi „(g, t). 



In order to prove this result, we will give a symmetric functions interpretation of our q, t- 
Narayana numbers: 



Theorem |6.1| , For all m > 1 and n > 1 we have 

Nara mjn (g,t) = (<^) m+n_1 • (Ve m+n _ 2 , /t m _i/t n _i), 

where e^ and hk are the elementary and the homogeneous symmetric functions of degree k 
respectively, V is the well known nabla operator introduced by Bergeron and Garsia (see [1, 
Section 9.6]), and the scalar product is the usual Hall inner product on symmetric functions. 

This result establishes a remarkable link between the q, i-Narayana polynomials and the well- 
known diagonal harmonics DH n , since Ve n is the Frobenius characteristic of this important 
module of the symmetric group & n , as shown by Haiman in [6] . 

Haglund [4J gave a combinatorial interpretation of the polynomial (V6 jn _|_ ri _2, h m ^ \h n — i) in 
terms of parking functions. In fact Haglund's result would be an easy consequence of the 
famous shuffle conjecture, which predicts a combinatorial interpretation of Ve n in terms of 
parking functions (see [51 Chapter 6]), if a proof of it could be found. 

In order to prove Theorem 6.1, we use the results of Section [5j proving that the combinato- 



rial polynomials in Haglund's result and our q, i-Narayana polynomials both satisfy the same 
recursion. This paper is organized in the following way: 

• In Section [2] we define three statistics on parallelogram polyominoes and two q, t- 
analogues of Narayana numbers. 

• In Section [3] we establish a bijection between our parallelogram polyominoes and a set 
of Dyck paths. We classify those words that are area words of members of Polyo mn . 
Area words are important in the definition of several statistics mentioned in Section |2] 

• In Section|4]we present a bijection from Polyo m n to Polyo n m which sends the bi-statistic 



(area, bounce) to the bi-statistic (dinv, area), thereby establishing Theorem 4.1 
In Section [5] we prove a recursion satisfied by both of our q, t-Narayana polynomials, 
which gives another proof of Theorem |4.1| 



In Section [6] we provide the necessary background to state Theorem 6.1, and we show 



how Theorem |6.2| follows from it. Theorem 6.1 is then proven. 



2. Three statistics on parallelogram polyominoes 

We may give an alternative characterisation of parallelogram polyominoes in terms of non- 
intersecting paths in the plane. This alternative charaterisation will prove useful in defining the 
statistics and mappings used in the paper. 

Consider a rectangular grid in Z 2 of width m and height n. On this grid consider two paths, 
both starting from the Southwest corner and arriving at the Northeast corner, travelling on the 
grid, performing only North or East steps, with the further restriction that they touch each 
other only at the starting point and at the ending point. Such a pair of paths uniquely defines 
a parallelogram polyomino. The region between the two paths is called the interior of the 
(parallelogram) polyomino. The two paths defining the parallelogram polyomino of Figure [I] 
are coloured in red and green, and the interior has been shadowed. 

In what follows we will encode a parallelogram polyomino as an area word consisting of 
natural numbers (unbarred numbers) and natural numbers with a bar on top (barred numbers), 
in the following way. We will label every North step of the upper (red) path with a barred 
number, and every East step of the lower (green) path with an unbarred number. This is done 
in two stages. 
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Figure 1. A parallelogram polyomino having a 12 times 7 bounding box. 



First, for each East step of the lower path we draw a line starting with the East endpoint and 
going Northwest until reaching the upper path: we label this step with the number of squares 
crossed by this line. Second, we label each North step of the upper path with the number of 
squares in the interior of the polyomino to the East of it which were not crossed by any of the 
lines that we drew during the previous stage. An example of this labelling is shown in Figure 
[2j where we put a black dot in the non-crossed squares. 
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Figure 2. The parallelogram polyomino of Figure [T] with its perimeter labelled. 

Once we have done this labelling, we read the labels in the following order: starting from 
Southwest and going to Northeast imagine moving a straight line of slope —1 over the polyomino. 
When we encounter vertical steps of the upper path or horizontal steps of the lower path we write 
the corresponding labels. If we encounter both types of steps at the same time then we write the 
label of the upper path first. The area word of the example in Figure[2]is 0112232221T211T2222. 

Notice that the sum of these numbers (disregarding the bars) gives the area of the polyomino, 
which is the number of squares between the two paths. This is the first of the statistics that 
are relevant to us. In the example the area is 30. 

Next we will define the dinv statistic. Consider the total order on the labels 

0<1<T<2<2<3<3<4<4<---. 

Given a polyomino with area word a\a2 ■ ■ ■ Ofc, we define its dinv as the number of pairs ctj, a,j 
with i < j and Uj is the immediate successor of a% in the fixed order. In the example of Figure 
[2J the number of such pairs containing is 4 and the dinv of the polyomino is 35. 

The last statistic that we introduce is the bounce. Consider the following path in a given 
polyomino: begin with a single East step from the Southwest corner, and then move North until 
reaching the East endpoint of a horizontal step of the upper path; at this point we "bounce" , 
i.e. we start moving East, until we reach the North endpoint of a vertical step of the lower 
path; at this point we "bounce" again, start moving North, and we repeat this procedure until 
we reach the Northeast corner. This path is called the bounce path. 

Once we have the bounce path, starting from Southwest corner, we label each step of the 
first sequence of vertical steps with 1, then each step of the second of such sequences with 2, 

3 



and so on; we label each step of the first sequence of horizontal steps with 0, then each step of 
the second of such sequences with 1, and so on. See Figure for an example of this labelling. 
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Figure 3. The labelled bounce path. 



The bounce of a polyomino is the sum of the labels of its bounce path, disregarding the bars. 
The bounce of the parallelogram polyomino in Figure [3] is 41. 

These three statistics give rise to a pair of bi-statistics on Polyo m n whose generating functions 

Nara m , n M):= £ 

PePolyo mn 

and 

Nara m , n (Q,t) := £ ^ounce(P) ^area(P) 
PePolyo mjn 

are studied in this paper. The polynomials Nara m , n (q,t) where first introduced in [3] by two 
of the authors of the present work. In the same paper, it was conjectured that these were 
polynomials symmetric in q and t, and as expressions symmetric in m and n. 

3. A BIJECTION WITH DYCK PATHS 

In this section we present a bijection ptd between the set Polyo m n and a set of Dyck paths 
having length 2(m + n). We then prove a result which shows how to read the area word of a 
parallelogram polyomino from its corresponding Dyck path under ptd. From this we will get a 
characterization of the area words of polyominoes from Polyo m n which will be used in the proof 



of Theorem 4.1 We finally observe that this description provides a way to computationally 
work with the set of area words of Polyo m n by working with the easier to construct set of Dyck 
paths. 

Recall that a Dyck path can be thought of as a path consisting of Northeast or Southeast 
steps lying between parallel horizontal lines, such that the path starts with a Northeast step, it 
never crosses the starting horizontal line, and returns to it at the end. Its length is simply the 
number of its steps it contains. Figure [4] shows an example of a Dyck path having length 38. 

Notice that a Dyck path is uniquely determined by the sequence of rises and falls we encounter 
as we move along the path from left to right. 

We will next describe a bijection between the polyominos in Polyo m n and the set of Dyck 
path of length 2(m + n) with m rises in even positions and n rises in odd positions, which do not 
return to the starting horizontal line until the end. This bijection appears in [2] in a somewhat 
different language. 

The idea is to read the steps of the upper and lower paths of a parallogram polyomino P 
alternatingly and form the Dyck path D = ptd(P) by using two rules. We perform a rise of 
the Dyck path for either a North step of the upper path or an East step of the lower path, and 
perform a descent of the Dyck path for either an East step of the upper path or a North step 
of the lower path. Using this construction, the polyomino in Figure [I] is sent to the Dyck path 
shown in Figure [4} 
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Figure 4. The Dyck path corresponding to the polyomino of Figure [TJ 

It should be clear that this mapping sends the parallelogram polyomines to the stated subset 
of Dyck paths. The fact that the Dyck path does not return to the starting horizontal line 
before the end corresponds to the fact that the upper and the lower paths do not intersect 
each other between the starting and ending points. The inverse operation is straightforward to 
describe and verify. 

We can easily read the area word of a parallelogram polyomino P from the corresponding 
Dyck path ptd(P) as we will now describe. Consider the Dyck path in Figure [4] when reading 
the next proposition. It consists of Northeast and Southeast steps lying between parallel lines 
which determine certain rows. 

Proposition 3.1. Let P £ Polyo mn with D = ptd(P). If we label the rows of D with 
0, 1, 1, 2, 2, 3, 3, • • • from bottom to top, then reading the labels of the rows of the rises from 
left to right we get the area word of the polyomino P. 

To prove this proposition, we will use induction on the number of pairs of steps of the upper 
and lower paths starting form the Southwest corner. At each step of this induction we will 
consider the partial box that includes the partial paths, i.e. the smallest rectangle that includes 
them (see Figure [5]). Then we imagine to complete the paths inside the partial box by moving 
along the edges to reach the Northeast corner, and we read the labels of the resulting polyomino 
on the partial paths. 



Figure 5 . The highlighted area is the partial box after the first 6 steps of both 
the green path and the red path. 

We claim that this gives exactly the corresponding part of the area word of the original 
polyomino. 

The key observation is the following claim. 

Claim. In the last pair of steps, the label of a North step of the upper path is always the distance 
from the right edge of the previous partial box, with a bar on top; while the label of a East step 
of the lower path is always the distance from the upper edge of the partial box. 

After proving this claim, it remains only to observe that the distance from the right edge of 
the previous partial box of the North steps of the upper path corresponds to the number of odd 
rows from the bottom line in the corresponding Dyck path; while the distance from the upper 
edge of the partial box of the East steps of the lower path corresponds to the number of even 
rows from the bottom line. This completes the proof of the proposition. 
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Proof of the Claim. At the beginning the upper path is forced to go North and the lower path 
is forced to go East. The partial box at this point consists of a single square, and we clearly 
have the partial area word 01: this is always the beginning of an area word for a polyonimo, 
and it corresponds to the first two rises in the corresponding Dyck path, as it should be. 

Now suppose that everything works up to a certain pair of steps, and let us make the next 
pair of steps. We have four cases (see Figure [6]): 




Case 1 



Case 2 




y y 
Case 3 Case 4 

Figure 6. The four cases. The previous partial box is violet, while the new 
one is black. 

Case 1 : The upper path moves East, and the lower path moves North. Then there are 
no labels to add, and the previous labels remain unchanged, since the partial box remains 
unchanged. 

Case 2 : Both the upper and the lower paths move North. Then the label of the North step 
of the upper path is clearly the distance from the right edge of the partial box, which is the 
same distance from the one of the previous partial box. The previous labels clearly remain 
unchanged. 

Case 3 : Both the upper and the lower paths move East. Then the label of the East step of 
the lower path is the distance to the upper edge of the partial box. The previous labels of the 
upper path remain unchanged, since in each row we added just a box crossed by the diagonal 
corresponding to the new East step of the lower path. The previous labels of the lower path 
also remain unchanged, since we did not move the upper edge of the partial box. 

Case 4 : The upper path moves North, and the lower path moves East. Then the label of 
the North step of the upper path is the distance from the right edge of the partial box minus 
1, since the first box becomes crossed by the diagonal of the East step of the lower path. But 
this is equal to the distance from the right edge of the previous partial box. The label of the 
East step of the lower path is clearly the distance from it to the upper edge of the partial box. 
The previous labels of the upper path remain unchanged, since in each row we added just a box 
crossed by the diagonal corresponding to the new East step of the lower path. The previous 
labels of the lower path also remain unchanged, since the diagonals of the previous horizontal 
steps all hit the upper path in the same spots as before. □ 

As an immediate consequence, we get a characterization of the words in the ordered alphabet 
0<1<1<2<2<3<3<-- - which are area words of elements of Polyo m n . 
We state this characterization here as a corollary. 
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Corollary 3.2. Consider the alphabet 0<1<1<2<2<3<3<---, with the letters in the 
given order. A word a±a2 ■ ■ ■ a r in this alphabet is the area word of an element of Polyo m n if 
and only if the following conditions hold: 

(1) a\ = 0, and this is the only that appears in the word; 

(2) there are exactly m of the a,i 's which are from the set of numbers without a bar {1, 2,3,...}, 
and exactly n of the a, 's which are from the set of numbers with a bar {0, 1,2,...} (in 
particular r = m + n); 

(3) for all i = 1, 2, . . . , m + n — 1, the letter a$+i is less than or equal to the immediate 
successor of the letter en, in the given order on the alphabet. 

We mention here that this bijection also gives an easy way to construct the polyomino from 
its area word: draw the corresponding Dyck path (this is immediate), and then look at the odd 
and even steps to construct the polyomino. 

4. The bi-statistics (area, bounce) and (dinv,area) 
This section is dedicated to proving the following theorem. 
Theorem 4.1. For all m > 1 and n > 1, 

Nara mj „(g,t) = Nara n , m (g, t). 

In order to prove this theorem it suffices to give a bijection from Polyo m n to Polyo n m which 
sends the bi-statistic (area, bounce) to the bi-statistic (dinv, area). 

The bijection that we will now describe is similar in spirit to the one used in the proof of the 
analogous (5j Theorem 3.15]. 

Let P 6 Polyo mn . Starting from P, we read the labels of its bounce path, getting a word 
consisting of barred and unbarred numbers. Then, starting from the bottom-left corner, for 
each turn of the bounce path, we look at the part of the path (upper or lower) that includes 
it. For example in the polyomino of Figure [3J the first turn of the bounce path is between 
and the next 1 in the labelling of the bounce path. The containing path consists of the first 4 
steps (counted from the Southwest corner) of the upper path. We label the vertical steps of the 
containing path with the labels used for the vertical steps in that part of the bounce path, and 
the horizontal steps of the containing path with the labels used for the horizontal steps in that 
part of the bounce path. See Figure [7] for an example. 






























































































1 






















1 






















1 



























Figure 7. The containing path and the new labels are blue. 

We then read the new labels by following the containing path from Northeast down to South- 
west. In the example we read 0111. 

During the remainder of the construction we will preserve the relative positions of these 
labels. 

We then repeat the algorithm with the second turn of the bounce path of P. In the example 
this occurs between the last 1 and the first 1 in the bounce path. This time the containing 
path consists of the steps of the lower path between the second and the eighth. We repeat the 
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procedure that we used before, and the word that we get reading the new labels will prescribe the 
relative positions of the l's and the l's. In the example (see Figure ^ we get the prescriptions 
1111111. This together with the other prescription gives a partial word 01111111. 

In general we will construct this partial word in a way that it can be the word of a parallel- 
ogram polyomino while respecting all the prescriptions. This will always be possible since the 
first step of the containing path that we read will always be labelled by the smallest of the two 
types of labels that we are considering: this is due to the definition of the bounce path. 
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Figure 8. The containing path and the new labels are violet. 



We keep doing this until all the labels of the bounce path of P have been included. At 
the end we will get a word of another parallelogram polyomino. In the example, at the next 
step we get the prescriptions 11211, which gives the partial word 011112111; then we get the 
prescriptions 222, which gives the partial word 01111222111; then we get the prescriptions 223, 
which gives the partial word 011112223111; then we get the prescriptions 3333, which gives 
the partial word 011112223333111; then we get the prescriptions 33443, which gives the partial 
word 01111222333443111; and finally we get the prescriptions 4444, which gives the final word 
01 1 TT222333444431TT. 



It is clear from the construction and the characterization of Corollary 3.2, that in this way 
we get the area word of a polyomino F{P) in Polyo nm . Moreover F{P) clearly has area equal 
to the bounce of the original polyomino P, again by construction. Figure [9] illustrates F(P) for 
when P is the polyomino of Figure [T| 
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Figure 9. The outcome of applying F to the polyomino of Figure [TJ 
We need to show that the dinv of F{P) is equal to the area of P. 

To see this, recall how we constructed the word of the new polyomino: for consecutive types 
of labels, we prescribed the relative positions by reading the corresponding containing path. 



But in the containing path, those pairs of vertical and horizontal steps which contribute to the 
dinv of the polyomino correspond each to a square in its area. 

It remains to show that F is a bijection. To see this, we can consider the inverse function: 
given a parallelogram polyomino, write in weakly increasing order its area word, and draw it 
as a bounce path with labels. Then reading the relative positions of consecutive types of labels 
you can reconstruct piecewise both the upper and lower paths. This completes the proof. 

Let us observe some remarkable consequences of this result. First of all, notice that iterating 
this bijection a second time, we get a bijection F o F from Polyo m n to itself which sends bounce 
to dinv. Moreover, applying the inverse and composing it with the flip along the Southwest to 
Northeast line that pass through the Southwest corner (which obviously preserves the area) we 
get a bijection from Polyo m n to itself which sends dinv to area. 

In conclusion, we see that all our three statistics are equidistributed both inside the same m 
times n rectangle and with the polyominoes in the flipped n times m rectangle. 



5. Recursions for Nara min (q, t) and Nara njm (g, t) 

In this section we prove that both Nara min (g, t) and Nara njm (g, t) satisfy a certain recursion. 
As an immediate byproduct we get another proof of the identity Nara min (g i , t) = Nara njm (g, t) 



stated in Theorem 14.11 

- — - — (r,s) 

Let Polyo m n be the set of polyominoes in Polyo m n whose labelled bounce path has r many 
l's and s many l's. In other words, r is the number of steps between the first and the second 
bounce of the bounce path, while s is the number of steps between the second and the third 
bounce. Define 



Nara 



ht) :- 



£bounce(P)^area(P) 



— — - (r,s) 

PGPolyo m „ 



(r,s). 



so that Nara mj „((7, t) is the sum over all r and s of Nara m n (q, t). Also, we define the (/-analogue 
of the non-negative integers by setting [0] q := 1, and for all positive integers n, 



n 



1 



l+q + q 2 + --- + q n ~ l . 



We define the q- analogue of the factorial of a non- negative integer by setting [0] q l := 1, and for 
all positive integers n, 



IB 



Finally, for < k < n, 



q [n-k] q \[k] q \ 
denotes the g-analogue of the binomial (^) . 

Theorem 5.1. For all m,n,r and s such that 1 < r < n and < s < m — 1, we have the 
recursion 

n—r m—s—l 



Nara^V) = f"+- 



h=l k=0 



s + r — 1 




's + h-1 


s 


q 


h 



M {H ' k) I +\ 



with initial conditions 



(n,s) 

Nara m,n (9.*) 



{qt)m+ n-X [m +n-2^ lfs = m _ l 

if s < m — 1 , 



(r,0) 



and Nara ln (q,t) = for r < n. 




IS 



m 



Figure 10. The black path is the bounce path. 



Proof. The argument in this proof is best understood by referring to Figure 10 

- — (r,s) 

The figure shows a typical element of Polyo, m . The orange grid cuts out an element of 

Polyo m _ s n _ r : its lower-left corner is placed at the beginning of the rightmost step of the 
bounce path labelled by 1. 

Observe that the labels of the bounce path in the orange grid are the same as the labels of 
the corresponding small path all increased by 1. Hence, together with the l's and the l's of 
the bounce path outside of the orange grid, we see that the bounce of the larger polyomino is 
m + n — 1 more than the bounce of the small polyomino in the orange grid. This shift is taken 
care of by the factor t m+n ~ 1 . 

The area of the larger polyomino is equal to the area of the small polyomino in the orange 
grid plus the yellow area, which is taken care of by the factor q r+s , the light blue area, which 
is counted by the factor f^ -1 ]^ and the pink area, which is counted by the factor [ s+ ^ ~ l ] ■ 
This explains the recursion formula. □ 

Let us denote by Polyo^m the set of parallelogram polyominoes in an n x m rectangle whose 
area word has r many l's and s many l's. Define 

)( g ,t):= Yl t area(P V inv(P) , 



Nara<$f 



PePolyo^ 



so that Nara rajm (q, t) is the sum over all r and s of Nara^j }(q, t). 



, ■ ( r ! S ) / 

These polynomials satisfy the same recursion satisfied by the Nara m n [q, t) s: 

Theorem 5.2. For all m, n, r and s with 1 < r < n and < s < m — 1, we have the recursion 

n-rm-s-l r in r , in 

.1 -I- r — I .1 -I- h. — I ii, i.\ 

S (<?,*), 



Nara£#( 9 ,t)=r* 



h=l k=0 



s + r — 1 




's + h-1 


s 


g 


h 



Nara 



(h,k) 



n—r,m—s\ 



(Qt) 



m+n—l rm+n-2] 
L m-l J 



if s = m 
if s < m - 



with initial conditions 

Narafe)(g,£) 

and Nara„2 (g, t) = for r < n. 

(r s) — 

Proof. Given an element of Polyon.'m with h many 2's and k many 2's, we construct an element 
of Polyo^l^ m _ s by subtracting 1 from all the letters in the area word, then removing all the 
resulting 0's and 0's and replacing the only —1 (which comes from the only 0) by 0. 

10 



For example, if we start with the word 0111221222111, which is an element of Polyo, 



(3,4) 
6,7 



with 



3 many 2's and 2 many 2's, then we first get —1000110111000, and hence we finally get 011111, 
which is an element of PolyOg^ 7 _ 4 = Polyog 3 ^. 

Now the area of this new element is clearly m + n — 1 less than the area of the original 
polyomino, since we subtracted 1 from all the letters of the area word different from 0. This is 
taken care of by the factor £ m + n_1 . 

The dinv of the original polyomino is equal to the dinv of this smaller polyomino, plus the 
dinv coming from the original and the l's, which is taken care of by the factor q r , the dinv 
coming from the l's and the l's, which is counted by the factor q s (the l's and the l's 

form a word which always starts with 1), and the dinv coming from the l's and the 2's, which is 
counted by the factor (as before, the l's and the 2's form a word which always starts 

with 1, but the dinv coming from this first letter is already counted by the that we insert in 
the new area word!). 

This explains the recursion. □ 

As already mentioned, these recursions give immediately Nara^m(g, t) = Nara mn (q, t), and 
hence another proof of the identity Nara mjn (g, t) = Nara„ )m (g, t). 



In this section we will use some tools from the theory of Macdonald polynomials. For a quick 
survey of what we need (and more), we refer to the book [1], in particular Chapters 3 and 9. In 
what follows we will recall only some basic facts, mostly to fix the notation. 

Let A = (J) n>0 A n be the space of symmetric functions with coefficients in C(q,t), where 
q and t are variables, with its natural decomposition in components of homogeneous degree. 
Recall the fundamental bases of symmetric functions: elementary {e^}^, homogeneous {h^}^, 
power {Pfj}^,, monomial {m^}^ and Schur {s^}^, where the indices [J, are partitions. 

A scalar product is defined on A by declaring the Schur basis to be orthonormal: 



where x ls the indicator function, which is 1 when its argument is true, and otherwise. Another 
fundamental basis of A is {H^}^, the modified Macdonald polynomial basis. 

The fundamental ingredient of the theory is the nabla operator V acting on A. This is an 
homogeneous invertible operator introduced by Bergeron and Garsia in the study of the diagonal 
harmonics DH n of & n . In fact, it turns out that Ve n gives precisely the bigraded Frobenius 
characteristic of DH n . 

The so-called shuffle conjecture predicts a combinatorial interpretation of Ve n in terms of 
parking functions. Special cases of this conjecture have been proven by several authors. In 
particular, Haglund [I] proved the combinatorial interpretation of (Ve n , hjh n -j) for 1 < j < n 
predicted by the shuffle conjecture. 

Surprisingly, this same polynomial provides the symmetric functions interpretation of our 
q, i-Narayana numbers. More precisely, we have the following theorem, which is the main result 
of this paper. 

Theorem 6.1. For m,n > 1 we have 



6. Symmetric functions interpretation 



> = X(A = n) 




Nara TO) „(g,t) = Nara n , m (g, t). 
li 



Moreover, we have 

Nara mj „(g,t) = Nara m , n (g,i) 



Proof of the Theorem 6.2, The symmetry in q and t comes from a general property of the nabla 
operator, which is an immediate consequence of the well-known identity [TJ Equation (9.8)]: 
nabla applied to any Schur function is symmetric in q and t. 



The second equation, symmetry in m and n, is obvious from the formula in Theorem 6.1 
Finally, the fact that Nara mj „(g, t) = Nara mjn (g, t) is a direct consequence of the symmetries 
and of Theorem 14. 1[ □ 



6.1. Proof of Theorem 6.1[ In order to prove Theorem 6.1 we need to make use of Haglund's 
combinatorial interpretation of (Ve m + n _2, h m -\h n -i). To do this we first require some defini- 
tions. 

For us a Dyck path of order k will be given by an area word which is a sequence of non- 
negative integers &1&2 '"bk such that b\ = 0, and b; L+ i < bj + 1 for all 1 < i < k. A domino is a 

pair of values (a, b) written as the first above the second 



A parking function PF of size k is a sequence of k dominoes 



ai 


a 2 




ak 


h 


b 2 







such that 



&1&2 • ■ • bk is the area word of a Dyck path of order k, and the a^s are a permutation of the 
integers {!,...,&} with the property aj < o$+i if frj < (and hence bi = bi + \ — 1). 



Example 6.1. PF 
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is a parking function of size 11. 



Remark 6.1. Parking functions are often represented by a diagram like the one in Figure 11 
In this diagram the red path represents the underlying Dyck path, where the number of the 
squares between the vertical steps of the Dyck path and the (green) diagonal are given by the 
lower numbers in the dominoes. The numbers that label the vertical steps of the Dyck path in 
the diagram are simply the upper numbers in the dominoes. 
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Figure 11. The parking function PF = 
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Given a parking function, we can reorder its dominoes by comparing first the bottom numbers, 
from the biggest to the smallest, and then, we place the dominoes with the same bottom number 
in order as we read them from right to left in the parking function. 

The reading word a(PF) associated to a parking function PF is the permutation that we 
obtain by reading the upper entries of this reordered sequence of dominoes. 

12 
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then we reorder the dominoes as 



Example 6.2. If PF 



and the corresponding reading word is a(PF) = 2 10 794811136 5. 

Given a parking function PF = ^ 2 °^ , we define its area to be area (Pi 7 ) 
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bi + . . . + bk, and its dinv d\nv(PF) as the number of pairs with 1 < i < j < k such that 
either bi = bj and a« < a,j, or bi = bj + 1 and aj > Oj. For example the area of the parking 



function of Example 6.2 is 14, while its dinv is 8. 

Given two disjoint sequences of numbers A and B, we denote by A LU B the set of shuffles of 
A and B, i.e. the sequences consisting of the numbers from A U B in which all the elements of 
A and B appear in their orig inal order, so that \A LU B\ = ( |A |^j B| ). 

For any a and b in N, we call Park a) b the set of parking functions PF of size a + b such that 
a(PF) £ (1, 2, . . . , a) LU (a + 1, a + 2, . . . , a + b). Finally, we set 



Parang, t):= £ ^-(PPyMPF) 

PFePark aj6 

We may state now the result of Haglund (see [1] for a proof, and [5] for the necessary back- 
ground) . 

Theorem 6.3 (Haglund). For all m > 1 and n > 1, we have 

(V 6m+n-2)^m-l^ii-l) — Pa^ra— l,m— 1 

This theorem reduces the problem of proving Theorem |6.1| to proving the following: 

Nara m , n (o, t) = (gir+^Para^i^i^, t). (6.1) 

In order to show the validity of this equation we do as follows. For 0<r<n, 0<s<m with 
r + s > 1, let Park^l^ m _ l be the set of parking functions PF of size m + n — 2 such that 

a(PF) e A LU B, \D (PF)nA\=r, and \D (PF) n B| = s, 

where ^4 = (1, 2, . . . , n — 1), B = (n, n + 1, . . . , m + n — 2), and Dq(PF) is the set of upper 
numbers of dominoes of P.F whose bottom numbers equal 0. Define the polynomial 



Para 



(r,s) 
n—l.m- 



^area(PF)^dinv(PF) 



PPePark 



(r,s) 

n — l,m — 1 



and set Para 



(0,0) 
n— X.m- 



_i(q,t) = 0. Clearly the sum of all these polynomials as r and s range over 
their possible values is equal to Para n _i im _i(q, t). With this new generalization in mind, it is 
clear that equation |6.1| holds true if 



(Qt) 



rn+n— 1 r 



Para^ m _ 1 (g,t) = Nara^)(g,t). 



Our proof of this identity is similar to what Haglund did in [4 J . 



We will show that (qt) 



m+n—l 



Para 



(r,s) 
n— l,77i- 



_i(q,f) also satisfies the recursion in Theorem 



with the 4-tuple (m, n, r, s) replaced with (n, m, s + 1, r), i.e. 



5.2 



(Qt) 



m+n— 1 



Para i- 5 l,m-l(9>*) 



m—s—1 n—r~l 

1 E E fl r 

fc=l h=0 



r + s 




'r + k-1 
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(gt) m+n - r - s - 2 Para 



(h,fc-l) (a t) 

n.-r-l,m-s-2\<i> l )i 



13 



which simplifies to 



Para 



(r,s) 
n— l.m- 



The initial conditions are 



Para 



(r,m— 1) 
n— l.m— 



iM) 



n— r— 1 m— s— 1 

£ £ 

/i=0 fc=l 



N(m,r) 
3f3m,n 



r + s 




r + k — 1 
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A; 



(h,k-l) 



(<?,*) 

(qt) m+n ~ l 




Para 



if r = n 
if r < n 



-X,m— s— 2 



($,*)• (6-2) 



and 



Para 



(r,s) 
0,m- 



Nara 



(s+1,0) 
m,l 



for s < m — 1. 



In order to see how the recursion |6.2| works for parking functions, we first make a sim- 
plification. It follows immediately from the definitions that, since the reading word of the 
parking functions we are interested in is a shuffle of the sequences A = (1,2,. . . ,n — 1) and 
B = (n, n + 1, ■ • • , n + m — 2), the pairs of dominoes with both upper numbers in A or both 
in B do not contribute to the dinv. The only pairs that contribute are the ones where one of 
the upper numbers is in A and the other is in B. Since all the elements of A are smaller than 
all the elements of B, we can simply consider dominoes in which the upper number is 1 (if the 
corresponding element was in A) or 2 (if the corresponding element was m B), with the dinv 
defined in the same way. 



For example the parking function PF £ Park 
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whose reading word is 



a(PF) = 16 14 7 8 15 4 10 11 5 12 6 13 1 2 9 3, 



would correspond to 
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whose reading word is 

a(PF) = 221121221212112 1. 
In both cases the dinv is 32, and the area is 18. 

Using this identification, what we do is the following: given an element PF of Park^_^ m _ \ , 
we remove the dominoes whose lower number is 0, and we decrease the lower number of the 
remaining dominoes by 1, keeping them in the given order. 

In our example, applying this procedure to PF we get 



r ( r > s ) 
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In doing this, observe that we will always get a parking function which starts with a domino Q 

which is always followed by a domino with lower entry 0, since PF cannot contain a sequence of 
three consecutive dominoes with strictly increasing lower numbers. This first domino contributes 
to both area and dinv and we can therefore remove it. In doing this we get an element of 



- rdIK n-r-l,m-s-2' 

dominoes in PF. 



where h is the number of 



dominoes in PF, and k is the number of 



Remark 6.2. Conversely, given an element of Paik n ^ r _ 1 m _ s _ 
domino, increase all the lower numbers by 1, and then insert r 
This gives us an element of Park^^ , 



2 , we can prepend it with a 



dominoes and s 



dominoes. 



_ x . In doing so, we are forced to put a 
14 



in front of the first 
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iii 



which we just prepended. Other inserted dominoes must satisfy the following constraints: a 

domino is followed by a domino 
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Let us now look at how the area and the dinv change with respect to the operation that we 
have just described. The area of the new parking function is equal to the area of PF minus 
(n — 1 — r) + (m — 1 — s), which is taken care of by the factor f n - r + m - s - 2 Q n the right hand 
side of 

The dinv is going to be the dinv of PF minus the dinv created by the dominoes that we have 

[2" 

removed. First of all, there are the pairs of dominoes Q 



creates dinv: this dinv is taken care of by the factor [ r 
there is the dinv created by the dominoes 1 
care by the factor [ 



and 



in PF, whose relative position 



on the right hand side of (6.2). Then 



in PF with the dominoes 



in PF: this is taken 



•+fc-ii 

k J 



since the first domino 



is necessarily preceded by a 



The initial conditions are obvious. This completes the proof of Theorem 6.1 
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